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Abstract 

We show nonlinear transport experiments on clean, suspended bilayer graphene that reveal a 
gap in the density of states. Looking at the evolution of the gap in magnetic fields of different 
orientation, we find that the groundstate is a spin-ordered phase. Of the three possible gapped 
groundstates that are predicted by theory for equal charge distribution between the layers, we can 
therefore exclude the quantum anomalous Hall phase, leaving the layer antiferromagnet and the 
quantum spin Hall phase as the only possible gapped groundstates for bilayer graphene. 
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The isolation of monolayer graphene [l[ has given a new twist to the research on two- 
dimensional electron systems, because graphene as a zero-gap semiconductor with a pseudo- 
relativistic dispersion relation shows fundamentally new effects like Klein-tunneling |2| that 
do not occur in conventional two-dimensional electron gases (2DEGs) with a finite bandgap 
and a parabolic dispersion. The presence of two atoms in the graphene unit cell is described 
by a sublattice-pseudospin, which acquires a Berry phase of 7r on closed k-space trajectories 
that include one of the K-points, which are the corner points of the first Brioullin zone and 
also constitute the Fermi surface of charge neutral graphene. One consequence of this Berry 
phase is an unusual Landau level spectrum that leads to the occurence of quantum Hall 
plateaus at conductances of half integer multiples of the Landau level degeneracy J], 4]. As 
the pseudospin degree of freedom is present also in Bernal stacked graphene multilayers, 
they can be described by one low energy theory that explicitly includes the chiral nature of 
the charge carriers [ij]. 

Bilayer graphene in A-B stacking is a 2DEG with a parabolic dispersion for small energies 
\E\ <C 71, where 71 ^0.4 eV is the interlayer hopping parameter that links two atoms sitting 
on top of each other 6| . As a consequence of the parabolic dispersion there is a finite density 
of states at the charge neutrality point (CNP) which, together with the weak dispersion leads 
to strong electron-electron interactions that make bilayer graphene at the CNP unstable 
towards interaction induced symmetry breaking 7- 1^. 

The exact nature of the electronic groundstate of bilayer graphene at charge neutrality is 
being discussed intensely at the moment, with numerous theoretical investigations suggesting 
a large number of phases. In each of these phases some of the three discrete degrees of 
freedom in bilayer graphene (spin, valley and layer) undergo a transition to a lower symmetry 
state. While it is clear that for large magnetic fields quantum Hall ferromagnetism will occur 



13|, and that for strong perpendicular electric fields a layer polarized state will form 3, Q], 



the nature of the groundstate at B=0 and Ei =0 is less obvious. Whereas some experiments 
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have found a conductive groundstate 
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161 ] . the majority of experimental results |l7H23| 



point to an insulating groundstate, an observation that limits the number of possible phases 



24| , bilayer graphene can be described 



to the ones that are bulk gapped. According to 
as a chiral 2DEG with chirality index J=2, the chiral degree of freedom being the sublattice 
pseudospin. 

In this communication we show experimental results that further elucidate the nature of 
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the bilayer graphene groundstate. We have performed nonlinear conductance measurements 
at the CNP of clean, suspended bilayer graphene samples as a function of magnetic fields 
oriented perpendicular and parallel to the graphene layer. Our measurements show a gap, 
that increases strongly only in a perpendicular magnetic field, whereas it stays constant to 
a high precision in a magnetic field that is exactly parallel to the bilayer plane. Our results 
allow us to make a statement about the spin order of the bilayer graphene groundstate 
and therefore to narrow down further the number of possible candidates for the electronic 
groundstate of bilayer graphene. 
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FIG. 1: a) SEM micrograph of the sample, b) Schematic of the sample structure, c) Conductance 
G as a function of backgate voltage V 9 . 



Bilayer graphene was exfoilated from natural graphite onto oxidized silicon wafers, and 
contacted with Cr/Au electrodes using e-beam lithography. The bilayer was suspended by 
partly removing the oxide in a buffered hydrofluoric acid wet etch and subsequent critical 
point drying, as described in 2l|. The data shown here were obtained on a four-terminal 
sample in Hall-cross geometry that had been shaped by an Ar/02 plasma dry etching process 
(figured]). We measured two-terminal conductance on contacts on opposite sides of the Hall- 
cross with the other two terminals floating or connected to a high-impedance voltmeter for 
Hall measurements. All measurements were done in a dip-stick dilution refrigerator at base 
temperature (T=60 niK). 

Figure [T] shows an SEM picture (a) and a schematic cross section (b) of the sample, 
together with a backgate characteristics of the device (c) measured at T=60 mK after 
current annealing 2l|. The charge neutrality point (CNP) is visible as a deep minimum 
in the conductance G measured as a function of gate voltage Y g at -0.6 V, indicating the 
formation of an insulating state around charge neutrality. The position of the CNP at a 



3 



small negative gate voltage indicates a very small amount of residual dopants on the flake. 
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FIG. 2: a) Conductance G as a function of perpendicular magnetic field and gate voltage V ff . 
b) cuts through a) at several fixed magnetic fields. 



Figure E^, shows a measurement of the zero-bias conductance G as a function of gate 
voltage V g and perpendicular magnetic field B^. From the line cuts in figure |2h one can 
see that there is a quantum Hall state at filling factor v—A developed for > LOT, giving 
a lower limit of 10 000 cm 2 /Vs for the charge carrier mobility. At higher fields additional 
intermediate plateaus start to occur, indicating a lifting of the spin and valley degeneracies. 

Differential conductance Grf=dI/dV at the CNP as a function of perpendicular magnetic 
field B± and bias voltage V s d is shown in figure [Hk- Several line-cuts at different magnetic 
fields are shown in figure |3b- Similar to previous studies [19l-l2l|. we find a strong suppression 
of conductance around zero bias, together with a BCS-like overshoot at a finite voltage of 
V s d ~ 3.5 mV, indicating the formation of an interaction induced broken symmetry state 
with a bulk gap. 
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FIG. 3: a) Differential conductance G^ at the CNP as a function of perpendicular magnetic field 
and bias voltage V s d- b) cuts through a) at several fixed magnetic fields. 



The Gd(Vsd) curves in figure Eb show two inflection points that we associate with two 
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different gap energy scales S and A 20|, |2l). Following the inffection points in finite perpen- 
dicular magnetic field, we find that the larger feature A increases linearly with 3.1 meV/tesla, 
whereas the smaller feature 5 shows a weaker response of 1.7 meV/tesla, both of them how- 
ever still much larger that the Zeeman energy for free electrons, which would amount to 
0.16 meV/tesla. 28] As we have shown in a previous study 20|, |2l) the larger feature A 
only exists in a narrow range of gate voltage around the CNP, whereas the small feature 5 
persists over the whole gate range that is accessible with the present sample. We think that 
5 is due to localization in disordered areas on the sample perimeter, and that A has to be 
associated with an interaction induced broken symmetry state located in the center region 
of the sample. 

Interaction effects in bilayer graphene have been described with the so-called broken 
symmetry state quasiparticle Hamiltionian 
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ley K and K', respectively, and A is the order parameter of the broken 
251. This Hamiltionian can be motivated by mean- field and renormal- 
.atioo g co Up calculation QflQ. factions w UI ,ead to a state wit h broken invert 
symmetry, with A pointing either along the z-direction, or perpendicular to e z in the x-y 
plane. In the former case a gapped state will result, whereas in the latter case a gap less, 
conducting phase with a structural anisotropy in the x-y plane, the so-called nematic phase 
occurs. In case of a gapped groundstate (A = (0,0, A 2 )) equation ([I]) predicts all electrons 
of the same spin and valley (the same spin-valley flavour) to be located in the same layer, 
a state that has also been described as a layer pseudospin magnet [7| • The order parameter 
A^, which is also called A in the theoretical literature, can depend on the electron spin 
s z and on the valley pseudospin t z , so that each of all four possible spin- valley flavours 
(t K, l K, I" K\\. K') can polarize to a different layer. This gives a total of 2 4 = 16 possi- 
ble gapped states, with 8 of them being qualitatively distinct and not related to any other 
by mere layer interchange. Of these 8 distinct states, 3 have no overall layer polarization, 
and can be assumed to be energetically favorable in case of vanishing external electric fields. 
These three states are the quantum anomalous Hall state (QAH), the layer antiferromagnetic 
state (LAF), and the quantum spin Hall insulator (QSH). Their symmetries are determined 



by different order parameters in the quasiparticle Hamiltionian (GQ), namely A z = Xt z for 
the QAH, A z = Xs z for the LAF, and A z = Xt z s z for the QSH. For the QAH, this has 
the consequence that electrons with different valley will sit in different layers, making the 
QAH a valley antiferromagnet, whereas LAF is a spin-antiferromagnet with the electrons 
in the bottom and top layer having opposite spin orientation. Due to a peculiar k-space 
topology 25[, each spin valley flavour has an intrinsic Hall conductivity, the direction of 
which is given by r z - sign(A 2 ): it depends on the valley and the sign of the order parameter 
in equation ([1]). While for the LAF and QSH states the (charge) Hall conductivities of the 
four different spin-valleys cancel out, in the QAH they add up, which means that for the 
QAH a nonzero quantized Hall conductance of 4e 2 /h is expected even in the absence of 
an external magnetic field. Although no net charge Hall conductivity is expected for the 
QSH phase, the resulting edge states will be helical, leading to a quantized spin Hall effect. 
More important, both QAH and QSH are expected to show a conductance of 4e 2 /h in a 
two-terminal transport experiment, in case the edge-states couple to the metallic contacts. 
However, as we have shown in a previous publication 21| . this condition is not necessarily 
fulfilled in experiment, as the clean part of a sample that hosts the gapped phase can be 
separated from the metallic contacts by a more disordered phase. This disordered phase 
will not let edge-states penetrate to the contacts, but will confine them to isolated puddles 
in the sample center. For the LAF charge, spin and valley Hall conductances cancel out, 
no edge-states can form, and consequently LAF is expected to be fully insulating even in a 
two-terminal conductance measurement. Note that the LAF state has also been predicted 
to show a gap whose size is independent of a perpendicular magnetic field, whereas QAH 
and QSH should show a strong increase of the gap by 8.5 meV/tesla [lj]]. The insensitivity 
of LAF to a perpendicular magnetic field has however been called into question recently, as 
this phase will seamlessly cross over to a quantum Hall ferromagnetic phase that shows a 
similar magnetic field dependence as QAH and QSH j^ . 

To identify the actual groundstate among the three candidates, it is useful to explore 
the three symmetries that might be preserved in them, namely time reversal symmetry 



T, spin rotation symmetry (SU2) and valley Ising symmetry (Z 2 ) 19|, l25j . In each of the 
three phases, only one symmetry is preserved: time-reversal symmetry T in the QSH, spin 
rotational symmetry SU2 in the QAH, and valley Ising symmetry Zi for the LAF. 

The valley Ising symmetry Z 2 , a state being invariant under exchange of K and K', is 
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broken in QAH and QSH, but is difficult to assess in a transport experiment. Time reversal 
symmetry T is broken by the QAH and LAF, and the case of QAH would manifest in a 
spontaneous Hall conductance at zero magnetic field, and in general as an inequivalence 
of G(B) and G(-B). In a real world sample however, puddles of different sign of the Hall 
conductance would exist next to each other, so that a spontaneous Hall conductance would 
be unobservable. Spin rotation symmetry SU2, which is broken in the QSH and LAF, with 
the spins in the bottom and top layer being locked into an antiferro magnetic arrangement, 
can easily be assessed by looking at the response of the gap to a parallel magnetic field, to 
exclude any orbital effects that are responsible for the big response to perpendicular field. 
While LAF is a real antiferromagnet, with the spins in the top and bottom layer pointing 
into opposite directions, the QSH also has to be considered a phase with antiferromagnetic 
spin order, however with the spin orientation being opposite for different valleys. In both 
LAF and QSH, the electron spins are not free to move, and a response to parallel magnetic 
field would be absent, as long as the Zeeman energy is smaller than the exchange interaction 
responsible for the antiferromagnetic spin arrangement j^]. In the QAH the electron spins 
are not ordered and would show a response to parallel field determined by the Zeeman 
spitting of the spin up and spin down levels. 
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FIG. 4: a) G^ at the CNP as a function of parallel magnetic field Bn and V s ^. b) Cuts through a) 
at several fixed magnetic fields. 



In order to assess the response of A to a parallel magnetic field, we have performed 
nonlinear conductance measurements in a vector magnet system that allowed to adjust the 
orientation of the magnetic field vector with high precision. As the Zeeman splitting of the 



conduction electrons with g=2 is much smaller than the response of A to a perpendicular 
magnetic field, the proper adjustment of the direction of By is crucial, and a small misalign- 
ment of a few degrees will lead to a perpendicular field component that will change A much 
more than what is expected from Zeeman splitting alone. Measurements were done with the 
sample mounted vertically, and Bii was applied in the direction of the bilayer plane. As ini- 
tial measurements showed an unreasonably large response of Gd(V s d) to By, we checked the 
correct alignment of the sample plane by applying an additional, small magnetic field in the 
horizontal direction, maximizing the height of the BCS-like overshoot visible in Gd{V s d) at 
V s d ~ 3.5 meV (see figures [3b and lib). We found that the sample had a slight misalignment 
of 1.6° with respect to the vertical direction, probably due to mechanical imperfections of 
the sample holder. Adding a small horizontal field component to the vertical magnetic field, 
we were able to compensate for this misalignment and to adjust the magnetic field to the 
sample plane with a precision of 0.1°. 

The results of a measurement of differential conductance at the CNP as a function of By 
and V s d are shown in figure HI Compared to the conductance as a function of perpendicular 
field (figure [3]) it is clear from figures Hh/b that the big gap A is not affected by a parallel 
magnetic field. Following the outer inflection points at V s ^ ~ 3 mV we find a magnetic 
field dependence of less than 70 /xeV/tesla, which is significantly smaller that the Zeeman 
splitting for an electron spin (160 /xeV/tesla). Our experimental results therefore allow us 
to exclude a change of A due to the Zeeman splitting of the conduction electrons. Within 
our experimental resolution, A is independent of a parallel magnetic field. Note that a 
small dependence of on B» occurs around zero bias V s d ~0. We think however that 
this change in Gd does not originate from the clean phase in the sample center and is not 
related to the gap A. According to the previous paragraph, this result shows that the 
electronic phase that we have labelled A is a spin-ordered phase, where the electron spins 
are bound by some exchange interaction and do not respond to an external magnetic field. 
Looking at the spin properties of the three possible candidates for A, we can conclude 
that A is not the quantum anomalous Hall (QAH) state, but that it has to be one of the 
remaining two: the quantum spin Hall (QSH) or the layer antiferromagnet (LAF). To further 
distinguish QSH from LAF one could make use of the spin Hall effect in the QSH, that 
should lead to a finite spin accumulation on the Hall terminals of a Hall cross sample. Such 
a spin accumulation might be detectable via the inverse spin Hall effect, using a metal with 



S 



strong spin-orbit scattering, like palladium, as contact material on the Hall terminals. Due 
to the non optimal homogeneity of multi-terminal, current annealed suspended graphene 



samples 21] . the significance of such an experiment would however have to be evaluated 
very carefully. 

In conclusion by performing nonlinear transport experiments on clean bilayer graphene for 
magnetic fields of perpendicular and parallel orientation we have shown that the groundstate 
at charge neutrality is gapped, and is a spin ordered phase. Of the three possible broken 
symmetry states that are suggested by theory, we can exclude that a quantum anomalous 
Hall state is realized in bilayer graphene. The groundstate has to be either the quantum 
spin Hall state, or the layer antiferromagnetic state. To distinguish the latter two from each 
other, further experimental work will be needed. 
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